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Abstract — The authors recently proposed a MIMO radar sys- 
tem that is implemented by a small wireless network. By applying 
compressive sensing (CS) at the receive nodes, the MIMO radar 
super-resolution can be achieved with far fewer observations 
than conventional approaches. This previous work considered the 
estimation of direction of arrival and Doppler. Since the targets 
are sparse in the angle-velocity space, target information can be 
extracted by solving an £± minimization problem. In this paper, 
the range information is exploited by introducing step frequency 
to MIMO radar with CS. The proposed approach is able to 
achieve high range resolution and also improve the ambiguous 
velocity. However, joint angle-Doppler-range estimation requires 
discretization of the angle-Doppler-range space which causes a 
sharp rise in the computational burden of the £± minimization 
problem. To maintain an acceptable complexity, a technique is 
proposed to successively estimate angle, Doppler and range in 
a decoupled fashion. The proposed approach can significantly 
reduce the complexity without sacrificing performance. 



I. Introduction 

Multiple-input multiple-output (MIMO) radar systems have 
received considerable attention in recent years. Unlike a 
phased-array radar, a MIMO radar [1] transmits multiple 
independent waveforms from its antennas. MIMO radar with 
widely separated antennas [2] exhibits spatial diversity, which 
improves target resolution. In colocated MIMO radar [3], 
phase differences induced by transmit and receive antennas 
can be exploited to form a long virtual array and thus achieve 
superior spatial resolution as compared to traditional radar 
systems. Compressive sensing (CS) is a recent development 
[4] -[5] and has already been applied successfully in diverse 
fields such as image processing and wireless communications. 
CS theory states that a if -sparse signal x of length N can 
be recovered exactly with high probability from O(KlogN) 
measurements via t\ -optimization. 

The application of CS to radar systems was investigated in 
[6]-[8], and to MIMO radar in [9]-[ll]. In [10], a uniform 
linear array was considered as a transmit and receive antenna 
configuration. In [9] and [11], the authors proposed a MIMO 
radar implemented by a small scale network. According to 
[11], spatially distributed network nodes, each equipped with a 
single antenna, serve as transmit and receive antenna elements. 
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The transmit nodes transmit periodic pulses. The receive nodes 
forward their compressive measurements to a fusion center. 
Exploiting the spareness of targets in the angle-Doppler space, 
an i\ -optimization problem is formulated and solved at the 
fusion center in order to extract target angle and Doppler 
information. This approach achieves the superior resolution 
of MIMO radar with far fewer samples than are required by 
conventional approaches. This implies low power consumption 
for the receive nodes. 

The contribution of this paper is a low complexity CS 
approach for obtaining range as well as direction of arrival 
(DOA) and Doppler information about the target. Range 
estimation cannot not be obtained with the method of [11], 
as the target range causes an identical phase shift to the 
signals received at all nodes and during all pulses. One way 
to obtain range information is to measure the travel time of 
the emitted radar signal [10]. However, the range resolution of 
such an approach might be limited. In this paper, we modify 
the transmitted waveforms of the scheme of [1 1], so that range 
information is reflected in different phase shifts over different 
pulses. In particular, the transmit nodes transmit pulses whose 
frequency increases by some step from one pulse to the next. 
In this way, we introduce the step frequency radar approach 
[12] to MIMO radar, an approach that results in high range 
resolution. Range resolution tends to increase as the transmit- 
ted signal bandwidth increases. However, wideband signals 
correspond to very short pulses, which experience low signal- 
to-noise ratios at the receiver. A step frequency radar transmits 
a sequence of longer pulses that are narrowband and together 
cover a wide frequency range. Thus, a step frequency radar 
transmits an effectively wideband signal that does not suffer 
from low SNR problems. An extension of [11] to joint angle- 
Doppler-range estimation would be straightforward; however, 
it would involve prohibitively high complexity. In this paper, 
we propose an approach to obtain angle-Doppler-range infor- 
mation in a decoupled fashion, which results in significant 
complexity reduction. In [13], a decoupled angle-Doppler 
estimation approach was proposed for the case of slowly 
moving targets. In that case, we can assume that the Doppler 
shift varies between pulses but within each pulse remains 
approximately constant. Based on compressively collected 
observations during one pulse, one can obtain initial estimates 
of the azimuth angles by discretizing the angle space only. 



Then, Doppler information is extracted by combining the data 
of multiple pulses. The basis matrix requires a discretization 
of the Doppler space only for the initial angle estimates. In 
this paper, we apply a similar idea to decouple the estimation 
of angle, Doppler and range. We propose to transmit a pulse 
train with constant carrier frequency, followed by a pulse 
train with carrier frequency that varies between pulses. Based 
on the received data during the first pulse train one can 
decouple angle and Doppler estimation along the lines of 
[13]. Based on these initial angle-Doppler estimates, the range 
information can be extracted from the data corresponding to 
pulses that have varying frequency. The proposed method 
significantly reduces the complexity as compared to the joint 
angle-Doppler-range estimation using CS without suffering 
performance degradation. 

II. Signal model for the constant carrier 

FREQUENCY 

Let us consider the same setting as in [11]. Assume K 
point targets and colocated antennas randomly distributed in 
a small area. The nodes transmit periodic pulses. The k-th 
target is at azimuth angle k and moves with constant radial 
speed v k . Let (r|, a*)/(r[, a\) denote the location of the i- 
th transmit/receive node in polar coordinates. The number of 
transmit nodes and receive nodes is denoted by M t and N r , 
respectively. Let dk(t) denote the range of the k-th target at 
time t. Under the far-field assumption, i.e., dk(t) ^> r-^ r , the 
distance between the zth transmit/receive node and the k-th 
target d\ k ld r ik can be approximated as 

d% r (t) * d k (t) - vt /r (O k ) = 4(0) - r£ /r (9 k ) - v k t (1) 



where r/- /r (6>/e) = r] /r cos(6 k - a/ r ). 

Assuming that there are Nj jammers located at (dj, 0j) and 
there is no clutter, the compressive samples collected by the 
l-th antenna during the ra-th pulse are given by 
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where 

1) c and / denote the light speed and carrier frequency, 
respectively; T is the radar pulse repetition interval; 

_ ■ 47Td k (0)f 

2) 7^ = /3 k e 3 c ; f3 k denotes the reflection coefficient 
of the k-th target; 

3) pimk = Vl(0k)f/c + f k {m - 1)T; f k = ^ is the 
doppler shift induced by the k-th target; 

4) IT S , I = 0, . . . , L— 1 represents the time within the pulse 
(fast time) and thus the pulse duration is LT S \ 

5) The i-th column of X contains the transmit waveforms 
of the length L from the i-th transmit node, where 
X^X = 1m ± \ 

6) <&i = 3hX^ (M x L) is the measurement matrix for 
the l-th receive node [11]; &i is an M x M t (M < M t ) 
zero-mean Gaussian random matrix; 



7) v{0 k ) = [ e ^^\...,e j ^^t^] T and D(/ fc ) = 
diag{[e^ 0T % . . . , eJ2*ML-i)T a ]y 9 

8) x jm denotes the waveform emitted by the j-th jammer 
during the m-th pulse; e/ m is the thermal noise at the 
l-th receive node corresponding to the m-th pulse; and 
the corresponding powers are fx- and a 2 . 

It can be easily seen that the phase term associated with the 

_ ■ 47rd fc (0) 

range e 3 a is independent of receivers and pulses and 
thus it can be absorbed by the reflection coefficient. 

III. Introduction of step frequency to MIMO 

RADAR USING CS 

Let us consider a MIMO radar system in which each 
transmit node transmits pulses, each of pulse repetition interval 
T, so that the carrier frequency of the m-th pulse equals 



/ m = /(l + A/ m ) 



(3) 



where A/ m is the frequency step, with 0<A/ m <l,m = 

N p . 

The baseband samples collected by the l-th antenna during 
the m-th pulse are given by 
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In ©, the phase term associated with d k (0) varies with the 
pulse index. 

Let us discretize the angle-velocity-range space on a fine 
grid: a = [(ai, b±, ci), . . . , (ajv, ^iv, cjy)]- Then © can be 
rewritten as 

N 9h f 

rim = V Sn e^-^X^D(^™)Xv m (a n ) + y Zm 



where s. 



_ / Ik, i 
n "l 0, c 



if the k-th target is at (a n , b n ,c n ) 



qir 



otherwise 

2c n / m Vi( a n)fm . 2b n f m (m - 1)T 



(6) 



and 
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In a compact matrix form we have ?j m = ^^ m s + y^ m , where 

Vim = ^X if [e^- 1 D(26 1 / m /c)Xv m (a 1 ), . . . , 

e^^D(2^/ m /c)Xv ro (a^)]. (8) 

If there are only a small number of targets as compared to 
TV, the positions of targets are sparse in the angle-velocity- 
range space, i.e., s is a sparse vector. A fusion center can 
combine the compressively sampled signals due to N p pulses 
obtained at N r receive nodes as 

f = [ff 1 ,...,ff JVp ,...,f^ JVp ] T = 0s + Y (9) 

where = [(* n ) T , • • • , (*i^ P ) T , ■ • • , (*n„n p ) t ] t and 
Y = [(yii) T >---.(yiJVp) r .---.(yJv r Af P ) T ] T - The vector s 



can be recovered by applying the Dantzig selector [14] to 
©. The location of the non-zero elements of s provides 
information on target angles, velocity and range. 



A. Unambiguous range and velocity 

In this section, we discuss the effects of step frequency 
on the unambiguous range R u and unambiguous velocity 
V u . In the case of slowly moving targets, i.e., j m ^T s L « 
1, k = 1, . . . , K, m = 1, . . . , N p , the Doppler shift change 
over the pulse duration T S L is negligible as compared to the 
change between pulses. Consider two grid points (di,bi,Ci) 
and (a J ? b^ Cj) in the angle- velocity-range space. Given = 
bj and Ci ^ Cj, there is no range ambiguity if 



, N p . It holds that 



[12]; 
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. If A/ m = 0, theni? M = f ; 
. If Af m = (m - 1) A/, then i? M - 
• If Af m is randomly generated within a predetermined 

range [/ m i n , fmax], then i? w -> oo when m is large. 
Similarly, let ai = a j9 6^ ^ bj and q = Cj. If 



e j4irbifm(rn-l)T/c _ e j4irbjf m (rn-l)T/c m 

velocity ambiguity will arise. Thus 

• If A/™ 

• If A/„ 



1, . . . , TVp, then 



0, then V u = 



(m — 1)A/, then equals the minimum 



, rn 



common multiple of { ^ T , 
If A/ m is randomly chosen from [/, 
14 — > oo when m is large. 



1,...JV P }; 



ai fmax\-> then 



B. Velocity resolution 

Next we investigate the effects of step frequency on the 
velocity resolution in terms of the column correlation in the 
sensing matrix. To simplify the analysis, we consider only one 
receive node. The sensing matrix for the l-th receive antenna 
is Si = ^J 2 . . . , ^in p ] T , where is defined in ®. 

On letting g& denote the z-th column of Si, the correlation 
of columns g& and g&/ equals 



Pkk' 



Em P =ivS(a*)Bf*v ro (a fc ) 
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fc = fe' 
(10) 



For simplicity, we make the following assumptions 

• To highlight the velocity resolution, let a/e = ay and 

ca; = Cfc/ ; 

• We consider the correlation of columns corresponding 
to the adjacent grid points in the velocity dimension, 
i.e., by — bj, = Ab. This is the maximum correlation 
of columns in the velocity domain, and thus dominates 
the velocity resolution. In this case, for slowly moving 
targets, B^' is approximately independent of 2bk J m ; 

• Assume Af m <C 1 and rj, i = 1, . . . , M t are sufficiently 
small. Then v m K) is approximately identical across 
pulses. 



Let a = 47rAbT f . Then dTUb can be approximated as 
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where A/ = [A/i, . . . , A/^J. 

A set of sufficient conditions that guarantee a reduction 
in the correlation of columns corresponding to adjacent grid 
points in the velocity dimension, are the following (see proof 
in Appendix I): 

A/jv p + i-m — m-1 ' L p/ J 

(note that A/ m = mAf satisfies this condition); 
• sin(an) > for n = 1, . . . , N p — 1. 
For example, for N p = b the sufficient conditions require 
that > | and a < \. It can be easily seen that a larger 
A/ m — A/Ar p +i- m , m > [iVp/2j can reduce the correlation of 
columns corresponding to adjacent grid points on the velocity 
axis. On the other hand, as A/ m ,ra = increase 
the bandwidth consumption increases. Therefore, there is the 
tradeoff between the velocity resolution and bandwidth. When 
Af m <C 1, the gain in the velocity resolution due to the 
introduction of step frequency is negligible. 

IV. Decoupled estimation of angle, velocity and 

RANGE USING CS 
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Fig. 1. Illustration of the decoupled angle-velocity -range estimation approach 

Solving an i\ minimization problem requires polynomial 
time in the dimension of s. Let us discretize the angle, 
Doppler and range space into [a±, . . . , a/vj, [&i, . . . , 6/v b ] 
and [ci, . . . , CjvJ, respectively. The joint estimation of angle, 
velocity and range requires complexity of O^NaN^Nc) 3 ). 
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For large N a , N b and N C9 the complexity cost of the CS 
approach would be prohibitive. In this section, we propose a 
decoupled angle- velocity-range estimation approach to lower 
the complexity burden. Let us consider the case of slowly 
moving targets, i.e., fkT s L <C 1. In this case, the Doppler shift 
within a pulse can be ignored. We propose to first transmit N c 
pulses with constant carrier frequency, and then transmit N s 
pluses with random step frequency. Based on these pulses, the 
estimation proceeds in the following steps. 

Step 1: For the data collected during one pulse, the phase 
terms associated with range and Doppler are constant over all 
receive nodes. Therefore, we can estimate the azimuth angle 
by discretization of the angle space only, as illustrated in the 
top graph of Fig. [T] and described in [13] in detail. The same 
process can be repeated on a number of different received 
pulses. Let I\ denote the set of angle estimates obtained 
based on the i-th pulse. The union of the angle-estimate sets 
Ti, . . . , Tn p are the angle estimates that are provided to the 
next step. 

Step 2: As described in Section [III the range information 
can be excluded from the basis matrix for the data with 
constant carrier frequency. Therefore, we can extract the 
Doppler information by applying CS to the data collected 
during the first N c pulses. The corresponding basis matrix can 
be formed based on the above obtained initial angle estimates 
and a discretization of the Doppler space as shown in the 
middle graph of Fig. [T] 

Step 3: By processing the received signal of N s pulses 
with the random stepped frequencies, we can extract range 
information as described in Section (TTTJ with the important 
difference that only the range axis needs to be discretized and 
used along with the angle- velocity estimates produced in Step 
2 (see the bottom graph of Fig. [T]). 

The complexity of Step 1 , Step 2 and Step 3 are respectively 
0(A(N a ) 3 ), 0((BN b ) 3 ) and 0((C7V C ) 3 ), where A, B and 
C are scalars much smaller than N a , N b and N c . Therefore, 
the total complexity of the decoupled scheme is 0(A(N a ) 3 + 
(BN b ) 3 + (CN C ) 3 ). For large N a , N b and N c , it holds that 
G(A(N a ) 3 + (BN b ) 3 + (CN C ) 3 ) « 0((N a N b N c ) 3 ) which 
implies significant savings. 

V. Simulations 

We consider a MIMO radar system with the transmit/receive 
nodes uniformly distributed on a disk of radius 10m. The 
center carrier frequency is / = bGHz and the pulse repetition 
interval is T = ^s. Each transmit node uses orthogonal 
QPSK waveforms. The received signal is corrupted by zero 
mean Gaussian noise. The SNR is set to dB. The SNR here 
is defined as the ratio of power of transmit waveform to that of 
thermal noise at a receive node. A jammer is located at angle 
7° and transmits an unknown Gaussian random waveform with 
amplitude 60. The target reflection coefficients are all one. 

Fig. compares the performance of velocity estimation with 
constant carrier frequency and randomly stepped frequency. 
The target scenes shown in Fig. [2] are generated via 100 
independent and random runs. The grey scale represents the 
times a target detected by CS occupies a particular grid point in 



the target scene. A lighter color indicates a higher occurrence 
frequency of a target. To highlight the velocity estimation, we 
consider an extreme case in which three targets are moving 
in the same direction of 0° and have the same range at the 
initial time, i.e. R = 1500m. The radial velocities of the 
three targets are 170m/ s, 17bm/s and 180ra/s, respectively. 
The unambiguous velocity for the constant carrier frequency 
/ = 5GHz is 120m/ s. The possible ambiguous estimates of 
three targets are bOm/s, bbm/s and 60m/ s, respectively. The 
number of transmit nodes and receive nodes are M t = 30 and 
N r = 5, respectively. M = 30 measurements are obtained at 
each receive node. The top graph shows the true target scene. 
We consider the worst case for velocity estimation in which the 
three targets are located at adjacent grid points in the velocity 
domain. Therefore, a single bright spot appears in the target 
scene instead of three spots. The last three graphs demonstrate 
the target scenes produced by the CS approach with constant 
and step carrier frequency. It can be seen from Fig. [2] that 
the introduction of random step frequency can eliminate the 
velocity ambiguity by using sufficiently many pulses (N p = 
10). In sharp contrast, the use of constant carrier frequency 
always yields the ambiguous estimates around bOm/s using 
the same number of pulses. 

Fig. [3] shows the estimates of target locations in the angle- 
velocity-range space using the proposed decoupled method. 
Three targets are moving in the directions of {—0.5°, 0°, 0.5°}. 
The radial velocity of the three targets are 70m/ s, 7bm/s and 
80m/ s, respectively. The corresponding ranges are 1200m/s, 
1250m/s and 1200m/s. We sample the angle- velocity-range 
space by the increment (0.5°, 5m/ s, 50m). The number of 
transmit nodes and receive nodes are M t = 30 and N r = 30, 
respectively. M = 30 measurements are obtained at each 
receive node. The transmitters first send N c = 5 pulses 
with constant carrier frequency and then N s = 5 pulses 
with randomly stepped frequency. The top and bottom graphs 
show respectively the true target locations and the estimates 
produced by the decoupled angle-Doppler-range scheme in 
100 random and independent runs. We can see that the 
information on the three targets is exactly recovered in each 
of the 100 runs. A false target arises in only one out of 100 
runs. The false target is quite close to the first target, i.e., 
(—0.5°, 80m/s, 1100m). The decoupled scheme requires only 
0.002% of the complexity of joint estimation of angle, velocity 
and range. 

VI. Conclusion 

We have proposed a CS based MIMO radar approach 
for obtaining angle- velocity-range estimates. First, we have 
introduced a step frequency approach in each transmit node, 
which not only achieves high range resolution but also improve 
the unambiguous velocity. Further, a decoupled angle- velocity- 
range estimation scheme has been proposed to alleviate the 
complexity burden of CS applied to the joint angle- velocity- 
range estimation. The proposed scheme can dramatically re- 
duce the computational cost of CS and still achieve good 
performance. 
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Appendix 

For a given pair (fc, fc'),fc ^ k f , satisfying a& = ay and 
H' — H = A6, h(Af) is proportional to the ratio of puw 
to pkk, which reveals the effect of Af on the correlation 
of columns corresponding to the adjacent grid points in the 
velocity dimension. Instead of analyzing ft, (A/), we define 
another function for convenience as follows: 



C(Af) can be expanded by the Taylor series of the first 
order as 



A/=0 



m=l 

N p Np + l-m 



= C(0)-2a E (A/ m (m-l) 

ra=LA/" p /2j+l n=m-l 

Af Np+1 _ m (N p - m)) sin(an) (13) 

Em=i ^fm d Q^^ |a/=o is required to be negative for 
C(A/) < C(0),A/ ^ 0. Therefore, sufficient conditions 
that guarantees C(Af) < C(0), Af ^ are 



> to m > [N p /2\ 



Af Np + 1 

(A/ m = mAf satisfies this condition); 



sin(an) > for n = 1, . . . , N p — 1. 
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Fig. 2. Velocity estimates using constant carrier frequency and randomly 
stepped frequency for M = Mt = 30 and N r = 5. 
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Fig. 3. Target scenes produced by decoupled angle- velocity -range estimation 
with randomly stepped frequency for M = M t = N r = 30 and N c = N s = 
5. 
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